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Quasi-bound states in periodically driven scattering 
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We present an approach for obtaining eigenfunctions of periodically driven time-dependent Hamil¬ 
tonians. Assuming an approximate scale separation between two spatial regions where different 
potentials dominate, we derive an explicit expansion for scattering problems with mixed cylindrical 
and spherical symmetry, by matching wavefunctions of a periodic linear drive in the exterior region 
to solutions of an arbitrary interior potential expanded in spherical waves. Using this method we 
study quasi-bound states of a square-well potential in three dimensions subject to an axial driving 
force. In the nonperturbative regime we show how eigenfunctions develop an asymptotic dressing 
of different partial waves, accompanied by large periodic oscillations in the angular momentum and 
a nonmonotonous dependence of the decay rate on the drive strength. We extend these results to 
the strong driving regime near a resonant intersection of the quasi-energy surfaces of two bound 
states of different symmetry. Our approach can be applied to general quantum scattering problems 
of particles subject to periodic fields. 


I. INTRODUCTION 


The main object of this paper is a time-dependent 
Schrodinger equation of the general form 


i<i> {r, t) = 




( 1 ) 


where the potential Vr is 7r-periodic in time (in units 
such that the fundamental angular frequency is 2, and 
h = m = 1). 

Various Floquet approaches have been developed for 
studying problems with similar formal equations in differ¬ 
ent parameter regimes, and Sec. [TTl gives a brief overview 
of relevant works. Here we present an approach for 
obtaining approximate eigensolutions of Eq. m , start¬ 
ing with explicitly known solutions for each separate 
Schrodinger equation with one potential Vj, j G {L,R}. 
Assuming a scale separation between two regions where 
either one of the two potentials dominates, quasi-periodic 
Floquet eigenfunctions of Eq. m are obtained by a time- 
dependent matching of the wavefunctions, described in 

Sec. mra 

As the main result of the paper, in Sec. lIIIBi we present 
a new explicit expansion for problems with mixed cylin¬ 
drical and spherical symmetry, 


i<j) (f, t) = 




Vnt(|r1)-F(i)-f 


0 {r, t) , (2) 


where Vnt is important up to some characteristic dis¬ 
tance from the origin, and F (t) is the periodic driving 
force. Relevant examples for such a geometry include, as 
detailed below, an atom in a linearly polarized laser field, 
a quantum particle interacting with a periodically-driven 
semiclassical scatterer (e.g. a cold atom interacting with 
an ion in a Paul trap), and interacting cold atoms or 
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molecules that are subject to oscillating fields. The gen¬ 
eralization to other settings, e.g. with a time-dependent 
potential in the interior region, or a potential with spher¬ 
ical symmetry in the exterior region, is straightforward. 
The presented matching conditions give immediately the 
complete spatial information on the wavefunctions. The 
explicit usage of analytic wavefunctions in each region 
gives access to fine details of the spectrum which may 
be hard to locate otherwise, and in particular the widely 
used Quantum Defect Theory (QDT) nm can naturally 
be used in the interior region. Our approach is nonper¬ 
turbative in both potentials Vj, but neglects the effect 
of either potential in some region of space. Therefore 
the obtained solutions can be considered, if necessary, 
as a starting point for an expansion that will treat the 
neglected contributions. 

Finally, in Sec. [IV]we employ our expansion to calcu¬ 
late solutions of Eq. (j2]) with a spherical square-well po¬ 
tential, and demonstrate general phenomena in the non¬ 
perturbative regime, e.g. nonmonotonous parametric de¬ 
pendence of the decay rate out of the well, large periodic 
oscillations of observables, and the resonant intersection 
of the quasi-energy surfaces of two bound states of dif¬ 
ferent symmetry. 


II. OVERVIEW OF FLOQUET SCATTERING 


If we consider —F (t) to be a monochromatic electric 
field amplitude, and Vnt (|^|) as the Coulomb potential 
for an electron, Eq. m describes an atom in an AC field 
(the AC-Stark effect), written in the length-gauge within 
the dipole approximation. In [^, Q it is proved that with 
Eq. ([2]), the bound states of Vnt (under general assump¬ 
tions) become resonances with an imaginary part which 
depends as a power-law on the amplitude of the pertur¬ 
bation - indeed the proof is perturbative in the electric 
field amplitude. 

This result is in fact general - for a Hamiltonian with 
a continuous spectrum of scattering states, the bound 
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states will generally turn into resonances under the ef¬ 
fect of a periodic perturbation (Ef. The reason is that 
the periodic perturbation makes every bound state with 
energy (— |c(;|) resonant with unbound states from the 
continuum of positive energy states, under absorption of 
at least n quanta from the perturbing potential (whose 
frequency is 2), where 

(-|cj|) + 2n>0, (3) 

and 2n gives the exponent of the power-law dependence 
of the resonance width on the perturbation amplitude. 

To be contrasted with the above picture, almost all 
of the bound states of a time-independent Hamiltonian 
can strictly survive the addition of a periodic perturba¬ 
tion, if the unperturbed Hamiltonian has a discrete s^c- 
trum (of isolated eigenvalues of finite multiplicity Q), 
a pure point spectrum (of discrete eigenvalues [7H9|), or 
a bounded continuous spectrum (in which case the per¬ 
turbation must obey certain conditions [ 13 ). The dis¬ 
creteness or boundedness of the spectrum in these cases 
stabilizes the spectrum (typically except under some spe¬ 
cific resonances with the external field), for any strength 
of the perturbation, and in these cases the proofs are 
nonperturbative. 

Nonperturbative studies of eqs. m-m have received 
a lot of attention within the intense-laser literature ini, 
where Eq. © is designated as being in the Kramers- 
Hanneberger (KH) frame. A variety of approaches have 
been developed for tackling this problem, focusing on 
different physical questions and in different parameter 
regimes of laser frequency, intensity and polarization. In 
this language, Eq. m expresses the fact that the rate of 
n-photon ionization is proportional to the n-th power of 
the field intensity, a result derived already in the early 
days of the field within Keldysh theory • 

A verj^eneral Green’s function approach was devel¬ 
oped in lillii and solved for Coulomb scattering of elec¬ 
trons, by neglecting all time-dependent terms except the 
leading-order averaged term. Known as the KH approx¬ 
imation, this approach is suitable in the regime where 
the frequency and intensity of the oscillating field are 
much higher than the atomic potential (in atomic units). 
In this limit, the perturbative picture of ionization rate 
which increases with intensity breaks, and the signifi¬ 
cant distortion of the effective (“dressed”) potential seen 
by the electron leads to the remarkable phenomenon of 
stabilization of the atom against ionization. Ref. [HI 
gives a comprehensive review of the works related to 
this effect. The interest in the “KH atom” has been re¬ 
newed in recent years following ex per imental results [Ti¬ 
ll^ and theoretical investigations jlfll . Recent works 
have also revisited the systematic expansion of an effec¬ 
tive time-independent Hamiltonian in the high-frequency 
limit [HI, 1^, an d the effects related to the potential’s 
initial phase [^ . 

Eor lower fields and frequencies, a wealth of techniques 
have been applied in the field. Photon absorption or 
emission processes which couple different scattering chan¬ 


nels were treated by numerically integrating the close- 
coupled equations [H, IHI, and employing a QDT ap - 
proach to extrapolate scattering cross sections [^, 1^ . 
Various effects related to this rich problem have been 
analyzed in simpler settings [28143^. including the ap¬ 
pearance and annihilation of bound states in the dressed 
potential, avoided crossing of resonances and their be¬ 
haviour in phase-space, resonant coupling between inter¬ 
nal levels, and nonmonotonic ionization rates (as laser 
intensity is increased). Within R-matrix theory in the 
Eloquet setting, space is divided into two regions and 
the solutions (typically obtained numerically within each 
region) are connected at the boundary [HI , and new ex¬ 
tensions of this approach have been recently suggested 
[39| . Numerical integration techniques have evolved in 
complexity and sophistication [HI and continue to be 
improved, and in particular, there is interest in calculat¬ 
ing and directly probing the angular distribution of pho¬ 
toelectron spectra [^, |H|- AC Stark shifts of trapped 
atoms have recently been modeled and measured in [^ . 

In a more general setting, the formalism for treat¬ 
ing Hamiltonians periodic in time, using Eloquet theory, 
is well known j43l-l45|. The periodicity of the Hamilto¬ 
nian allows defining an extended Hilbert space in posi¬ 
tion and time, in which the scalar product is defined to 
include integration over the temporal period. This ex¬ 
tended Hilbert space can be spanned by set of spatially 
orthogonal wavefunctions and a Eourier basis for time- 
periodic functions, e.g. by all wavefunctions of a spe¬ 
cific TT-periodic Hamiltonian. Then, any other 7r-periodic 
Hamiltonian can be expanded using such a basis for the 
extended Hilbert space, and all of the tools of time- 
independent quantum theory are available, which can be 
powerful in many scenarios, e.g. for employing perturba¬ 
tion theory. Time-dependent perturbation theory is also 
widely used [dbl , and usually the interest is in transition 
rates between asymptotically time-independent states. 

To conclude this section, we briefly mention the re¬ 
cent interest in atomic systems in the ultracold regime 
(i.e. with temperature < 1 mK, 113 ), which are trapped 
by oscillating fields. Overlapping a trap for neutral atoms 
with a periodically driven Paul trap for ions [isl, was 
suggested in [49| and realized first in [Eol, followed by 
the demonstration of a trapped ion immersed in a di¬ 
lute atomic Bose-Einstein condensate [El|, Ell ^ many 

other experiments. The effect of the periodic drive of 
the ion has been analyzed for classical collisions with the 
atom [EEf, for quantum scattering employing a master 
equation description [EJ, and for an ion and atom in 
separate traps [HI- As mentioned before. Quantum De¬ 
fect Theory (QDT) constitutes one of the most common 
theoretical tool for modelling the short-range part of the 
interaction in atomic scatterring setups, and continues to 
evolve [EEMEEI. together with new models and methods 
|6EI63, applied to many-body states as well [HI, [HI • As 
a last example for a driven system with two-bocK inter¬ 
action we mention polar molecules in AC traps [TQI, UM • 



3 


III. QUASI-BOUND STATES IN and the normalization applicable to a square-integrable 

TIME-PERIODIC POTENTIALS wavefunction is 


In this section we formulate a method for finding wave- 
functions of Eq. o, with a time-independent potential 
Vl (f) which is assumed to be significant inside some in¬ 
terior region |r| < d, and a time-dependent 7r-periodic 
potential Vr (r, t) which dominates in the exterior region 
|r I > d. The essential assumptions at the basis of the pre¬ 
sented approach are that the wavefunctions of each of the 
potentials can be found explicitly, and that there is some 
meaning to dividing space into the interior and exterior 
regions, even if only as a (zeroth-order) approximation. 
We focus in this paper on finding quasi-bound states (res¬ 
onances), while scattering states follow the same expan¬ 
sion with just a redefinition of the unknown parameters. 

A. Matching conditions for Floquet-expanded 
wavefunctions 

To simplify the basic expressions, we take a one¬ 
dimensional (ID) notation for the derivation in this sub¬ 
section, starting with the ID equation corresponding to 

Eq. (H}, 


0* (x, t) 0 (x, t) dx = 1. (7) 

Expanding the functions at the matching point we write, 

4^L,uj-\-2n ~ ( 8 ) 

= (9) 

k 

SO the first matching condition of Eq. implies 

E «2nC2ne-'""* = E (10) 

n j,k 

which gives 

C2n<^2n — ^ ^ ^2(n—j)^2jf • 

3 

A similar expansion for the gradients, 

'^^'L,uj+2n (d) = / 2 n, ( 12 ) 


icj) (x, t) = 




VL{x) + VR{x,t) 


4>{x,t). (4) 


To formulate the matching conditions of the quasi- 
periodic wavefunction (j) (x, t) at a boundary point d 
which separates the regions of the left potential Vr {x) 
and the right potential Vr (x,t), we consider the ansatz 


f En a2n4>L,oj+2n (x) x <d 

4>ix,t) = < 

[ En ^■2n4>\u,+2n X >d 

(5) 

where 0 L,u;+ 2 n {x) and uj+ 2 n (1^^ latter being tt- 
periodic) are solutions of the Schrodinger equation with 
potential Vr (x) and Vr (x, t) respectively, and energy 
(quasi-energy) cj + 2n. Eor notational simplicity, the 
above summation does not indicate explicitly a summa¬ 
tion over any degeneracy of the wavefunctions, which 
must involve independent matching coefficients. In the 
ID case this can include left- and right-going waves (if 
the boundary conditions allow), and in higher dimensions 
there could be summation over other quantum numbers. 
It is also assumed here that some prescribed boundary 
conditions at x ^ ±oc are already included in (t)R^uj+ 2 n 
and uj+ 2 n' omit the explicit range of Eourier sum¬ 
mations on integers n G Z. The wavefunction (j) (x, t) 
is parametrized by uj which can be chosen in the range 
—2 < cj < 0, however uj does not determine uniquely the 
wavefunction - there can be different functions with the 
same value of uj (but different coefficients). 

The matching conditions at d are 


(jjR (d) = (l)R (d, t), V(/)L (d) = V(I)R (d, t ), (6) 


^rR,.,+2ni.d,t) = Y,92ke-^^’^\ (13) 

k 


gives 


f2nCi2n = '^g2(n-j)hj- (14) 

3 

The two matching relations can be written in matrix 
form (once a finite truncation has been applied), 

Cd = Db, Fa = Gh. (15) 

where a, b denote the expansion coefficients in vector no¬ 
tation and C^D^F^G are matrices. By writing the two 
equations in block form 

the compatibility of the two matching conditions implies 
the vanishing of (at least one) eigenvalue (or, more gener¬ 
ally, singular value in the SVD decomposition) of K (cj), 
with the corresponding kernel vector then giving the ex¬ 
pansion coefficients. The same arguments can be applied 
to the smaller matrix (since G and F are assumed diag¬ 
onal in the current expansion, and would in general be 
invertible) 

FC-^Db = Gb ^ {G-FC-'^D)b = 0, (17) 

whose kernel vectors give the exterior region coefficients 
6, from which a immediately follows. 
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Defining 


2k,2n — 


POO 

/ (^> 0)* <Pk^+2n 0) 

Jd 


and similarly, 


fd 

Q2k,2n = / dX(t)l^^+2k {X, 0 )* (t>l,u:+2n 0 ) > 

J — oo 


(18) 


(19) 


we get for the normalization condition (which can be eval¬ 
uated at t = 0), the bilinear expression 


/ oo pd poo 

(t)*<t)dx= / (t)*(t)dx+ / (t)*(t)dx 

-OO J — OO Jd 

— ^ ^ ^2kQ‘^k,2n^2n ^ ^ ^2/c~^2/c,2n^2n ~ I 5 (^0) 

k,n k,n 

or in matrix form, 

d'<Qa + PPb = l, (21) 

SO that normalization can be guaranteed by dividing a, b 
by the square-root of the l.h.s. We note that this nor¬ 
malization is relevant only if the entire wavefunction is 
square-integrable, and we will discuss the case of wave- 
functions with free-particle components in the next sec¬ 
tion [following Eq. 


B. 3D matching of wavefunctions with mixed 
cylindrical and spherical symmetry 


In this subsection we write the matching conditions for 
the 3D problem 


(r, t) = 


1 


V" + M„t(|r1)-F(t).f+V(t) 


(t> {r, t) 
( 22 ) 

with a spherically-symmetric potential in the interior re¬ 
gion, and a linear (periodic) drive in the exterior region, 
where Vi (t) is added for convenience as detailed below. 

For the general linearly-driven time-dependent 
Schrodinger equation (in the exterior region) 


# (f, t) = 




(p (f, t), (23) 


a family of solutions can be written in the form 
(/)(r,t) oc e^<r(t)-r-ig(t) ^ 

with 

q{t) = F + k 


(24) 


(25) 


where k is the (possibly complex) constant of integration, 
and 


9it) = + 


f 


k.Fit') + -F{tf+VAt') 


dt'. 
(26) 


Below we will need the cylindrical waves, i.e. the solu¬ 
tions of the time-independent free particle Hamiltonian 
in cylindrical coordinates (p, (p), defined by 

(f; k,a) = ^kpsma) (27) 

where m is the magnetic quantum number, k is the 
wavenumber which can in general be complex, and a is a 
complex parameter. is a Hankel function of the 

first or second kind (corresponding to outgoing and in¬ 
coming traveling waves respectively), or a Bessel function 
(which we denote with a superscript J). 

Specializing to the case that the drive is 7r-periodic and 
coaxial at any time, we can choose a fixed cylindrical 
coordinate system in which 

F{t) =F^{t)z. (28) 

We will further simplify the current expressions by taking 
Vi (t) of eqs. ([2^ - ([23|) to cancel the A?-independent term 
in Eq. ([26|) , so that outgoing and incoming traveling-wave 
solutions to Eq. (|23]) can be written using Eq. (|27|) in the 
form 


Jn’m {r,t-,k,a) 

^ ^-i^k‘^t^iF^{t)z^{^,2) j'- 1 ^ 29 ) 


Fixing the magnetic quantum number m which is con¬ 
served, and defining k 2 j = y^2 (cj -f- 2jf), the most general 
wavefunction that solves Eq. (|23|) for the exterior region 
|r| > d is 


(pR,m {r, t)= Fj j ^2i, O) , 

d J Ca 


a=l,2 


(30) 

which takes at each value of k 2 j a superposition of out¬ 
going and incoming cylindrical waves, parameterized by 
integrals in complex a-plane along two contours Ca with 
weight functions h^ 2 j both to be determined in the 
following. 

The above expansion becomes useful by using a rep¬ 
resentation of the spherical Hankel function of the first 
kind as an integral over cylindrical waves 0, [73| in the 
form 


hf-'^ {kr)Pp {cose) 

= [ dasma^{-iy~"' Pl^ {cos a) a) (31) 

JCi 

where P[^ are the associated Legendre polynomials and 
the directed contour of integration Ci depends on k. For 
k with a positive imaginary part we must take Ci = 7r/2+ 
i (oc, —oo), and (kr) then decays asymptotically as 
e~\^\r/r (we note that {f]k^a) diverges at p ^ 0, 
however the integral, which gives {kr), is well defined 
for any p > 0, and decays for r ^ oo, which is just what 
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we need). For k real and positive the contour of integra¬ 
tion is given hy Ci = i (oc, 0 ) + [ 0 , ir] {ir i ( 0 , —oo)}. 

This directs us to take the arbitrary weight function 
for outgoing waves in Eq. (|3Q|) to be of the form 

b 2 j («) = - (cosa), (32) 

h 

with 62 j,Zi become matching coefficients and the con¬ 
stants are determined by boundary conditions at 

infinity as detailed below. In App. [^we show that each 
term in the summation of Eq. (|32]) . when plugged into 
the integral over Ci [in Eq. (|3Q|) ]. can be written in the 
following form; 


/ da sin ab^ {a) (f, t; k 2 j, a) = 

JCi 


—e 




E (33) 


ll,l 


where are spherical harmonics, and the outgoing and 
incoming radial functions in the exterior region are de¬ 
fined by 


E ^h,h,h,h,ijh {F"' it) k2j)ju it)rj (k2jr) 

(34) 


with the coefficients cq^/ 2 ,/ 3 ,U,/ being defined in Eq. dsn. 


For any value of fc, we have similarly to Eq. m 


ji {kr)Pr (cose) = 

f dasina^ P[^ (cosa) , (35) 

J[o, 


where ji is a spherical Bessel function, and Eq. (|35]) allows 
to express the spherical Hankel function of the second 
kind (for k with nonnegative imaginary part) by using 

= 2ji — Therefore, using the fact that the 

/, m-dependent coefficients in Eq. m and Eq. (|35]) are 
identical, we can replace in Eq. (|3Q|) the integral over the 
contour C 2 (which need not be further specified) by an 
expression identical in form to the expansion in Eq. (EH), 
with the outgoing waves replaced by (minus) incoming 
waves. The wavefunction expansion in 3D analogous to 
Eq. ® is then 


with 


(t>m (r, t) 


'En,l a 2 n,ie-'(“+ 2 ”)VL, cD+2n,z {r ) Yj^ ( 6 >, ip) \r\ < d 




?( 2 ) 


( 1 ) 


?(1) 


(36) 


(37) 


In the interior region the wavefunction (j)L,uj+2n,i (^) is a solution with energy cj + 2n of the Schrodinger equation 
with potential Vl (F), expressed in spherical coordinates {r^O^cp), and the summation includes all partial waves and 
energies. In the exterior region we have incorporated the effect of boundary conditions at infinity into the 

choice of which will be detailed at the end of the current section. 

Therefore, equating the wavefunction on the surface of the sphere \f\ = d we have 


^e-*(“+2")*a2n,;<(>L,c.+2n,; id)Yr = E-^Wi^ id^t)Yr 

n,l j,li I 

= E E = E ^ 62 ,,qd 2 yq,/, 2 (n-,Wr (38) 

j,li l,p n,l j,li 

where are the expansion coefficients of the Eourier series of (t)\ 2 j,h,i which in general must be obtained 

numerically. Then the first matching condition is 


0L,u;+2n,Z {d) a2n,l 


Yd2j,li,l,2(n-j)i>2jM- 

j,h 


(39) 
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Since we equate the interior and exterior wavefunctions on the surface of a sphere, we need to consider only the radial 
derivatives. This gives 


dr(t)L,oj+2n,l (d) a2n,l = ^ ,/,2(n-i) , 

jM 


(40) 


with g 2 j,ii,i, 2 p the expansion coefficients of the Fourier series of drcj)^ 2 j h i {d^t). The latter derivative can be written 
explicitly as 


^rR 2 j,h,l (^’0 ~ (t) ^ 2 j) 

X [f- {t)ju+^ (f- {t) r) h[f ik 2 jr) +ju [f- (t) r) [i (h + h) h[f [k^jv) - k^jh^^Xi (^2i0] } • 


(41) 


For the normalization integrals performed at t = 0 as in Eq. (1201, we find using the orthogonality of , 

rd 

^ 7 -U \ ■ 1 \ / r 

^R, 2 j',lXK 


_^ pd _ _^ poo ^ 

1 = E “2n',Z«2n,i / drr^ 4>l^^+2n',l4>L,uj+2n,l + Y ^F2j,lul- (42) 

n,n' 0 ' I 


In App. [O we lay down for completeness the expan¬ 
sion of integrals which are required in order to calculate 
expectation values of general rank-0, -1, and -2 tensor op¬ 
erators (we restrict the expressions to axially symmetric 
wavefunction with m = 0). The expectation value of any 
time-independent (or 7r-periodic) operator is always tt- 
periodic for the Floquet eigenstates. The normalization 
integrals in Eq. (j42]) are a special case of Eq. (jC5p , which 
can be evaluated at t = 0 if the wavefunction is square- 
integrable. However, for expansions which contain free- 
particle components (discussed in the following), when 
we can only integrate over the bound components of the 
wavefunction, the normalization integral is 7r-periodic be¬ 
cause the relative weight of the nonnormalizable compo¬ 
nents oscillates in time. In this case one must divide 
expectation value integrals by the squared norm, both of 
which being 7r-periodic functions. 

In the interior region, the values of the above integrals 
can be obtained without explicitly performing the inte¬ 
gration, directly from the wavefunctions and their gra¬ 
dients at the matching point. This can useful especially 
when the interior wavefunctions are not explicitly known 
close to the origin, but rather are determined within a 
QDT formulation [Tj • We will use the following gen¬ 
eral notation for wavefunctions in spherical coordinates; 

cl>{r) = Yir)Yrie,ip), (43) 

r 

such that u (r) is a radial function obeying a one¬ 
dimensional Schrodinger equation with an effective po¬ 
tential which includes the centrifugal barrier. The pro¬ 
jection of two eigenfunctions and 02 of the interior 
Hamiltonian with energies ei and 62 correspondingly, is 
shown in App. [B]to be given by 

J = (^1 — ^ 2 )”^ ~'^ 2 '^i| I 5 

( 44 ) 


where = drUi^ and 5 i, 62 are assumed to have equal 
imaginary parts. The left-hand side of Eq. (|44|) gives the 
integrals required for the normalization, with the factor 
of 2 relevant for the off-diagonal projections (when 0 i 7 ^ 
02 )- In the limit of 0i ^ 02 we have 


f 


r^dr = - lim (ei - £ 2 )' 

Z Sl^S 2 


U1U2 — U2U1 


which gives the diagonal normalization terms. 


d ’ 

(45) 


Einally, we comment on the constants in 

Eq. (|3^ . which should be chosen in the following way, 
consistent with the domain of validity of Eq. (|3T]) . Eor 
the partial waves with fHe {cj} + 2 j < 0 the wavevector 
k 2 j can be chosen with a positive imaginary part and 

setting = 0, h ~ ^ gives exponentially de¬ 

caying outgoing waves. In the terms with 94e {ci;} + 2j > 0 
two types of boundary conditions at infinity can be im¬ 
posed. If we assume that positive energy states are free 
(and form a continuum), we should set S^ 2 jii ~ 

h ~ ^ search accordingly for resonances with 
Jmlcj} > 0 and 3m {k} > 0, which correspond to in¬ 
coming waves whose amplitude diverges at infinity. The 
imaginary part of uj gives the rate of formation of the 
resulting quasi-bound state. Otherwise, if there is some 
potential at infinity which reflects waves inwards, solu¬ 
tions will have uj purely real, we can set h ~ 

^R 2 jii the relative phase of waves reflecting from 

the boundary, assuming that it depends only on the en¬ 
ergy and the partial-wave angular momentum quantum 
number li of the nondriven problem. 
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IV. LINEAR DRIVE WITH A SQUARE-WELL 
INTERACTION 


0 

-1 


In this section we employ the methods presented in the 
previous sections to study a model system consisting of a 
spherically-symmetric square-well potential and a time- 
dependent periodic linear drive which acts outside of the 
well. We demonstrate the analysis of general phenomena 
in the nonperturbative regime. It is interesting to note 
that quantum wires and dots have been modeled by 
similar finite-barrier potentials, and the expansion pre¬ 
sented here can be used to solve a mixed-type system. 

Using the frequency of the periodic drive, Q, we can 
define the length and energy scales 



• s-wave 

— p-wave 

— d-wave 


9 




8.6 


do = y^2h/mQ^ Eo = hQl2, (46) 

and the variables become nondimensional by rescaling 
according to 

T — y vjdo^ k — y kdoi t — y tQ/2, 

after which we have explicitly h = m = 1 and the drive’s 
frequency in these units is Q = 2. With a spherical 
square-well potential, 

Kvell (|r|) = 1 0 °’ 

(where d and Vb are nondimensional, measured in the 
units of Eq. (|46]) ). the Schrodinger equation in the inte¬ 
rior region becomes 


i4) 


-^V' + Kvelldfl) 


(49) 


and the regular solution inside the well is a spherical 
Bessel function. 


FIG. 1. The spectrum of bound states in the square well 
potential [Eq. (liSl) ] with nondimensional width d = 1.15, as a 
function of the depth Vb- There are two s-waves, one p-wave 
and one d-wave. The least-bound (s-wave) state reaches the 
threshold at Vb ~ 8.5 (not shown). In figures [2] to [H we study 
the properties of this bound state when a periodic driving 
force with frequency U = 2 is turned on. 



(r, t) oc {kr) (50) 

In those units, we take the periodic force of Eq. (|28|) to 
be a simple harmonic drive with amplitude F 2 , 

(t) = F 2 cos 2t. (51) 

In Eig. [1] we show the spectrum of the time- 
independent square-well over a small range of Vb values 
at d = 1.15. Eor these parameters, the external drive 
fixed at frequency Q = 2 does not resonate with any of 
the transition frequencies between the states. The least 
bound s-wave state which we study in the following is 
pushed towards the threshold at Vb ~ 8.5. 

Eigure [ 2 ] shows the imaginary part of the Eloquet 
quasi-energy of the resonance state which evolves from 
the least-bound s-wave as a function of F 2 and Vb- hol¬ 
lowing the discussion at the end of the previous sec¬ 
tion, we solve for the state with positive imaginary part, 
that gives (half) the transition rate for formation of 


FIG. 2. Imaginary part of the Floquet quasi-energy of the 
state which evolves from the least-bound s-wave state of 
Fig. [U as a function of the periodic drive amplitude and the 
potential depth. A nonmonotonous dependence on the pa¬ 
rameters can be seen, where beyond some Vb-dependent value 
of F 2 , an increase of the drive amplitude leads to a decrease 
of the decay rate of the state out of the well. This decrease 
terminates when the imaginary part of uj reaches zero, which 
happens along a continuous line in parameter space (that here 
appears “jittered” due to the discrete sampling). Beyond this 
line the resonance seems to have disappeared from the quasi¬ 
energy spectrum. 


such quasi-bound states. Since the driving force is time- 
reversal invariant, taking the complex conjugate of this 
state and reversing the sign of t gives the solution which 
describes quasi-bound states decaying out of the well. We 
will therefore (somewhat loosely) refer to the imaginary 
part of the quasi-energy as the decay rate. 

For low drive amplitude, the quasi-bound state’s decay 














FIG. 3. Expectation value and standard deviation of the 
quantum average of the squared angular momentum opera¬ 
tor P (nondimensional, with h = 1), of the driven quasi¬ 
bound s-wave state. The parameters of the state correspond 
to the highest point in Fig. [ 2 l with Vo = 8.6 and F 2 = 0.235, 
and the state’s quasi-energy has shifted from uj ^ —0.0225 
(see Fig. [J) to cj —0.0612 + 0.00155z. The large amplitude 
oscillations are coherent, while the large standard deviation 
indicates that the state is an asymptotic broad distribution 
of different partial waves. 

rate grows parabolically (as can be inferred from a log- 
log plot, not shown) which is the expected perturbative 
result [Eq. (|3])]. In the nonperturbative regime the decay 
rate is clearly nonmonotonous; for a strong enough drive 
the decay rate begins to decrease and then reaches zero, 
where the two complex-conjugate resonances meet and 
seem to annihilate and be removed from the quasi-energy 
spectrum. We note that this is not a threshold effect - 
the real part of the quasi-energy is separated from 0. 

In Fig. [3] we show the expectation value and the stan¬ 
dard deviation of the quantum average of the squared an¬ 
gular momentum operator [calculated using eqs. (jC5p - 
(|C6|) ]. for the driven quasi-bound s-wave state at the 
highest point in the parameter region of Fig. [2l Both 
quantities display a large amplitude oscillation over one 
period of the drive - we note that this oscillation in it¬ 
self is coherent and involves no uncertainty. On the other 
hand, the probability distribution of the angular momen¬ 
tum (at any fixed time within the period) is seen to be 
very broad. We can infer that the fact that the expec¬ 
tation value remains close to zero is misleading (/ is not 
a good quantum number even approximately), and un¬ 
der the effect of the drive the state develops a strong 
superposition of many partial waves. We note that the 
imaginary part of the energy (which gives a decaying ex¬ 
ponential envelope) is ignored here. 

The nature of this superposition can be further seen 
in the solution coefficients of the expansion in Eq. (|3^ . 
which are depicted in Fig. 0] for the same state. The 
quasi-bound s-wave state which for F 2 = 0 would have its 
entire amplitude at (n = 0, / = 0) and {j = 0, /i = 0), has 
developed a broad superposition of partial waves (mostly 



FIG. 4. Absolute value of the matching coefficients of the 
solution ansatz of Eq. ()36D . in the interior region {a 2 n,i) and 
in the exterior region ( 52 j,q), for the state of Fig. O Outside 
of the range of (n,/) and (j,/i) values shown the coefficients 
quickly decay. We note that although inside the well the 
coefficients correspond exactly to a partial wave expansion, 
in the exterior region the parameter h does not correspond 
directly to a single partial /i-wave for F 2 / 0 . 



FIG. 5. The Floquet quasi-energy uj (mod 2 ) of two bound 
states in the square potential [Eq. (l48l) ] with nondimensional 
width d = 2.3 as a function of the depth Vq. The s-wave 
lies in fact deeper in the well than the p-wave for the entire 
parameter range - its real energy is a; — 2 , and the “crossing” 
of the states (calculated for F 2 = 0 ), becomes relevant only 
when considering F 2 7 ^ 0 , when the periodic drive resonates 
with the energy difference of the two states, as treated starting 
with Fig. [1 


outside of the well). The “checkboard” pattern is the re¬ 
sult of the dipolar nature of the coupling, which conserves 
(- 1 )”+' [or (-iy+'"]. 

Starting with Fig. [5] we consider a well which is twice 
wider and supports more bound states. The quasi¬ 
energies UJ (mod 2) of two of the bound states are plotted 
in this figure as function of Vq, around a point of cross¬ 
ing. The energy of the s-wave state is lower by —2, so 
that the crossing is only of quasi-energies (mod 2), and is 
irrelevant for the time-independent well. Figure [6] shows 
the real part of the quasi-energies for the same states, in 
dependence on both Vb and F 2 in a small region of pa¬ 
rameters. A resonance line emanates from the crossing 
point at F 2 = 0, on which the periodic drive’s frequency 
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FIG. 6. Real part of the Floquet quasi-energy a; (mod 2) of 
the two states of Fig.[5]in dependence on the depth of the well 
and the strength of the drive, around the crossing of quasi¬ 
energies. Emanating from the crossing point is a singular 
line of resonance (a “seem” of the two surfaces), on which 
the periodic force mixes completely the two states. The blue 
surface (lower on the left of the seem, upper on the right of 
it) corresponds to the state evolving from the s-wave bound 
state, and the red surface is the p-wave. 



FIG. 7. As in Fig. [6l only showing the imaginary part of the 
quasi-energy. The s-wave surface is blue, and remains below 
the p-wave surface for most of the parameter region except in 
a small part to the right of the seem. 


resonates with the energy difference of two states. This 
line constitutes a singular “seem” of the two quasi-energy 
surfaces (we note that the surfaces can be trivially con¬ 
tinued into the F 2 < 0 region). 

Figure [7] shows the imaginary part of the quasi-energies 
of those two states, which is also plotted on a log-log scale 
in Fig. [HI On the line of degeneracy the two states mix 
completely. Away from this line, the decay rate of the 
s-wave, whose real part of the energy lies in the range 
—4 < ^zuj < —2, grows quartically with F 2 , while the p- 
wave’s decay rate grows quadratically (and is also much 
larger), which is again the perturbation theory result. 
However, close enough to the resonance, it can be seen 



FIG. 8. As in Fig. (T] only showing the imaginary part of the 
quasi-energy on a log-log plot. The red surface of the p-wave 
is now upper for most of the parameter region, with an almost 
constant slope equal to 2, the perturbation theory result of 
quadratic decay rate out of the well for this state. The deeper 
bound s-wave has slope of 4 except within a Vb-dependent 
distance from the resonance, where its strong mixing with 
the p-wave lowers its power-low exponent to 2. 



FIG. 9. The quantum expectation value of the squared- 
angular momentum operator for the same states as in fig^ [5ll8l 
with the range of Vo somewhat larger. The value of (L^) is 
obtained by averaging over a period of the driving force. The 
red surface of the p-wave starts at / (/ + 1) = 2 for F 2 = 0 
and quickly decreases as a function of F 2 , depending on the 
distance to resonance. 


that the strong partial mixing of the s-wave with the 
p-wave changes the dependence of the former on F 2 to 
quadratic. The threshold of this regime passes already 
at F 2 < 0.01 for values at the edges of the figures, and 
goes down (towards F 2 0) closer to the resonance. 

Figures [9] and [To] show the expectation value of and 
its standard deviation, for the same s- and p-wave states 
over a somewhat larger range of Vb values around the 
resonance. For F 2 = 0 the value of the nondimensional 
angular momentum for the s-wave is 0, and for the p- 
wave it is / (/ + 1) = 2. The standard deviation of the 
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FIG. 10. As in Fig. [9l but showing the standard deviation of 
the squared-angular momentum operator (period-averaged). 
The p-wave (red) surface grows steeply to its maximal value of 
1, achieved after full mixing with the s-wave state, along two 
lines which emerge from the crossing point in both directions. 
Beyond these maxima lines, the p-wave is pushed further into 
an s-wave inside the well, while developing an asymptotic su¬ 
perposition of different partial waves. These two lines can 
also be identified in Fig. 1111 and can be considered to border 
the nonperturbative regime. 



FIG. 11. As in Fig. (TUI but showing the pe^-to-peak ampli¬ 
tude of the periodic oscillations of Stddev (T^), and the axes 
have been rotated for clarity. The two lines of null oscillation 
amplitude of the standard deviation correspond to the max¬ 
ima lines in Fig. [10] and dehne the nonperturbative regime. 


quantum distribution is 0 for both. 

The angular momentum expectation value of the p- 
wave is seen to quickly go down as it is mixed with the 
s-wave for F 2 7 ^ 0 , and close enough to the resonance, this 
happens for arbitrarily small F 2 . At the same time the 
standard deviation peaks at 1 for the p-wave (at some Vq- 
dependent value), the result of it being mixed with the s- 
wave and strongly pushed into the well. Then, for further 
increase of F 2 the p-wave turns more and more into s- 
wave within the well, while developing a superposition of 



FIG. 12. Absolute value of the matching coefficients of the 
solution (as in Fig. in, for the state continued from the p-wave 
of hgs. [MTTl at Vo = 9.1076 and F 2 = 0.0003. We note that 
for F2 = 0 this state has its entire amplitudes concentrated 
at (n = 1, / = 1) and {j = 1, /i = 1). Already at such a weak 
drive (but very close to resonance), we find (P) ~ 0.0004 and 
Stddev (L^) 0.36, as the state is just entering the nonper¬ 

turbative regime with amplitudes of new partial waves becom¬ 
ing important in the solution, in addition to the two resonat¬ 
ing bound states. With further increase of F2, the state will 
develop a notable admixture of partial waves, mostly outside 
of the well. 

The amplitude of the periodic oscillations in the stan¬ 
dard deviation of is shown in Fig.[TTJ With the driving 
force being relatively weak {F 2 ^ 0 . 02 ), these oscillations 
(in time) have a small peak-to-peak amplitude (compare 
with Fig. However, the peak-to-peak amplitude itself 
shows a nontrivial oscillatory pattern in its dependence 
on the parameters, which is the result of interference of 
partial waves with different values of j and therefore dif¬ 
ferent factors of . 

The null lines of the peak-to-peak amplitude, which 
emanate from the crossing point in both directions in 
the parameter space and coincide with the maxima lines 
of the mean standard deviation, can be taken as delimit¬ 
ing the nonperturbative parameter region. In Fig. [12] we 
show the matching coefficients for a state which is very 
close to resonance and at a very weak drive starts to en¬ 
ter the nonperturbative region. The amplitudes of new 
partial waves in addition to the two resonating bound 
states of the time-independent well begin to contribute 
to the solution, as the drive is strong enough to local¬ 
ize the higher energy p-wave inside the well and quickly 
dress it with an asymptotic tail of more partial waves. 
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Appendix A: Expansion of linearly-driven cylindrical waves in spherical waves 

The proof of Eq. (j33j) proceeds by using Eq. ([29]) to write 
^iF^{t)z da sin axm^ (E; k 2 j,a) {a) 

= dasinaxm {f;k2j,a) i‘2j,hPh (cosa) (-d^ ( 2 h + l)jh {-F"" (t) k2j) Ph {cosa) 

= (-*)'" ( 2^2 + 1 ) W ^ 2 i) E ,3 W , PiT) Ic, da sin (f; k 2 j,a) Pl^ (cos a) 

= b 2 j,iJh {F^ W k 2 j) El, {k 2 jr) (cos 0 ) 

= Eu ( 2^4 + l)ju {F it) r) Pu (cose) Ei„i, b 2 j,iJi 2 {F^ W ^ 2 ,) E ;3 {^ 2 ^) PJ^ {cosO) 

= EhmMM W ^ 2 i)i 4 (i)^) (fc 2 ir) E; ci,M„u,iYr iO,p) 

(Al) 

where the multiplicative factors have been omitted for simplicity, and by using the definition of 

I (^, t) given in Eq. (|M|) . Eq. (jAip results in Eq. (|33]) . In the derivation of Eq. (jAip . the plane-wave expansion in 
terms of spherical Bessel functions has been used (twice), the coefficients of expansion of a product of two (associated) 
Legendre polynomials (which can be written using Wigner 3-j symbols) are defined by 


W , P™=, P;7) = P™^ (w) P™^ («;) P;7 {w) dw, 

the coefficients cq^/ 2,/3 obtained using Eq. (|3T]) and Eq. (jA 2 p and given by 

= 2 { 2 I 2 + 1 ) (- 1 )'^ , P,^, P^) , 

and the coefficients Qi,Z 2 ,/ 3 ,U,z similarly 

CiMM = CIM ( 2^4 + l)i^^W {Pl^,Pu,Pr) /Nr, 

with the definitions 

r™ (61, = Nppr {cos 9) 


Nr = (-i)™E(ETT) 74 ^^(rE^Wi^TH!. 


(A 2 ) 

(A3) 

(A4) 

(AS) 


Appendix B: The projection of two eigenfunctions of the internal Hamiltonian 


In order to derive Eq. (|44|) . let € 1,62 be the (possibly complex) energies of two complex eigenfunctions 0i,02 of 
the interior Hamiltonian Hl- Eor the projection of the two within the interior region, we can write (since both sides 
vanish) 

0 = {(/2 \{Hl - £ 1 ) \h ) - (<(>1 \{Hl - £ 2 ) \h ). (Bl) 


By canceling the potential energy terms, we get after rearranging the kinetic terms and terminating the integration 
at an arbitrary point d (which is allowed since the equality above holds identically in space), 

{eicplcpi - e 2 (p*i 4 > 2 ) r'^dr = {u^d/ui - uld^U 2 ) dr 


1 

2 


{u^drUi — uldrU2)\^ — Jq {drU^drUi — drUldrU2) dv 


(B2) 


where the factor of 1/2 appears since we assume that the nondimensional kinetic energy term is — In the second 
line of the above equation, the integrated term is purely imaginary being the difference of two complex conjugates. 
Taking the complex conjugate of the entire equation and adding, this term drops and we get 


/■ 


2 91e{(5i — S 2 ) = —-2 91e {u 2 dr 


which gives immediately Eq. (lUl) . 


UldrU 2 }\a, 


(B3) 
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Appendix C: The expectation value of tensor operators 


In this appendix we give explicitly the expansion of integrals which are required in order to calculate expectation 
values of general tensor operators, in the Floquet eigensolutions of Sec. IIIIBl For simplicity we treat here only the 
most useful case of axially symmetric wavefunctions, with m = 0 (no ip dependence). 

Using the notation introduced in Eq. (|43j), we start by writing the 7r-periodic part of the wavefunction in the form 

r (r,t) = {r,t)YP (Cl) 

n,l 


which corresponds to the expansion in Eq. (j36j) of wavefunctions in the interior region. Eor such wavefunctions, we 
define the (unnormalized) expectation value in the interior region of a purely radial operator O (r), 




%[ 0 {r)] = I dVir (r,t)|"0(r) = 

1,1' n,n' 


^2n,Z^2n' 


J d^a [U2 




(C2) 


The above expression can be rewritten as 


3o[0{r)]=J2li,i[0{r)], (C3) 

I 

where we have defined for convenience the functional 

Ii,i' Pir)] = E J dr [u^ 2 n,iT^ 2 n'A , (C4) 

with the summation taken over pairs of states enumerated by {(n, /), (n', /')} with fixed I and I'. 

Eor example, the normalization integral calculated for any time (generalizing (j42j) for expansions with free-particle 
components) can be written using the above notation as 


3o [i]=E^M [i], 

I 


(C5) 


with i the identity operator. Any other expectation value must then be divided by the value of this normalization 
integral. Similarly, the expectation value of the squared angular momentum operator is given by 


3^0 




— E +1) di,i [i] 


(C6) 


For an operator of a general radial part multiplied by the position vector, O (r) r, only the Cartesian z-component 
survives the integral (for axially symmetric wavefunctions), and we can write using zjr = cos9 

3i [O (r) r ] = [ d^f\r (r, tp O (r) f = E^M' E [ dr *Oir)r (C7) 

1,1' n,n' 

with the coefficients being 

Pi^if = 2 ttNi Ni, J dO sin 0 cos 0Pi {cos 0) Pif {cos 0). (C8) 

Using that fact that nonzero terms will have |/ — /'| = 1, we find 


3i [O (r) f] = [O (r) r]. (C9) 

I 

Eor an operator with a general radial part multiplied by a bilinear combination of position vector components, 
O (r) where a, /3 G {x, y, z}, only the diagonal terms with a = f3 survive the integration (for wavefunctions with 
m = 0), with the result 


32 [O (r) fafp] = E 9a, I,I'll,I' [O (r) r^] , 

1 , 1 ' 


(CIO) 
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where 


qa,i,i’ = 27rA^i 


/ dO sin 0 


cos “1“ 2 ^ {,^a,x H” ^a,y^ 


Pi {cosO) Pi' {cosO). 


(Cll) 


In all of the above expressions, Ii^i' [O (r)] as defined in Eq. (jC4p is valid in the interior region. To get the complete 
result for expectation values in whole space, the integration over the exterior region must be added, where the 
wavefunctions are expanded differently in Eq. (|3^ . In this case, Eq. m is to be replaced by 


(r, t) = Yj Yj (r, t) Yi, 

j,h I 

and accordingly, Eq. (El becomes in the exterior region 
lix[0{r)]= Y 

n n' U V 


(C12) 


(C13) 
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